Probing Dark Energy through Scale Dependence 
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We consider the consequences of having no prior knowledge of the true dark energy model for the 
interpretation of cosmological observations. The magnitude of redshift-space distortions and weak- 
lensing shear is determined by the metric on the geodesies of which galaxies and light propagate. 
We show that, given precise enough observations, we can use these data to completely reconstruct 
the metric on our past lightcone and therefore to measure the scale- and time-dependence of the 
anisotropic stress and the evolution of the gravitational potentials in a model-independent manner. 
Since both dark matter and dark energy affect the visible sector only through the gravitational 
field they produce, they are inseparable without a model for dark energy: galaxy bias cannot be 
measured and therefore the distribution of dark matter determined; the peculiar velocity of dark 
matter can be identified with that of the galaxies only when the equivalence principle holds. Given 
these limitations, we show how one can nonetheless build tests for classes of dark energy models 
which depend on making measurements at multiple scales at a particular redshift. They are null 
tests on the model-independent observables, do not require modeling evolution in time and do not 
require any parametrization of the free functions of these models — such as the sound speed. We 
show how one can rule out or constrain the whole class of the most-general scalar-tensor theories 
even without assuming the quasi-static limit. 
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I. INTRODUCTION 

The introduction of dark energy relaxes many of the 
assumptions behind ACDM cosmology. Not only can the 
background evolution be altered, but in many models 
the dark energy is predicted to have significant pertur- 
bations. As a result, the simple ACDM relations between 
the dark matter perturbations and the gravitational fields 
sourced by them are no longer valid. The detection of a 
form of dark energy different from a cosmological con- 
stant hinges on the unambiguous determination of how 
these relations are broken. 

Various approaches to this goal have been taken. On 
the level of background evolution, null tests for the com- 
patibility of the cosmic expansion history with ACDM 
were proposed in Rcfs [1, 2], while a principal-component 
analysis of the equation of state constraints was de- 
veloped in Ref. [3]. On the level of linear pertur- 
bations, parametrizations are usually used to give a 
model-independent description of the growth of large- 
scale structure [4-8] , although principal component anal- 
ysis has also been performed [9] (see [10] for a short 
review). Alternatively one may try to focus on under- 
standing the limits that a general-relativity-like theory 
must obey, but describing the remaining freedom with 
a parametrization framework [11-13] or attempt to de- 
scribe the properties of dark energy through an effective 
framework [14-17]. Such approaches, when contrasted 
with data, usually require parametrizations in order to 
break down degeneracies, but which simultaneously in- 
troduce paramctrization-dependent biases. 

Our approach, introduced in Ref. [18], has been to take 



a step back and ask what the observables which we can 
extract from cosmological probes are when we do not as- 
sume a particular model for dark energy. We have argued 
that as a result of the presence of unknown dark energy 
perturbations, the only model-independent observables 
are ratios of quantities typically considered to be unam- 
biguous measurements. This implies that such variables 
as the dark-matter growth rate / or dark-matter per- 
turbation amplitude as are not actually observable in a 
model-independent way. 

In this paper, we elucidate further this framework, 
showing that a helpful way to understand our previ- 
ous conclusions is to realize that our cosmological probes 
(galaxy counts and weak lensing) rely on tracers which 
propagate on geodesies and therefore can at best deter- 
mine the gravitational potentials which affect their mo- 
tion. Indeed, we show that one can in principle map 
out the whole metric on our past lightcone and there- 
fore that the slip parameter describing anisotropic stress 
is a direct observable not requiring any parametrization. 
On the other hand, what we cannot do at all is to ever 
probe the dark matter distribution using these means, 
since the dark matter does not interact directly with our 
probes, only by sourcing the gravitational potential. This 
means the effective Newton's constant, which is impor- 
tant to quantify dark energy perturbations, cannot be 
observed. Given this difficulty, we nonetheless describe a 
strategy for constraining dark energy models, generaliz- 
ing our work away from the quasi-static limit considered 
previously. 

We start in section II by defining the notation, and 
reviewing the assumptions and the model-independent 
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observables that we have found before. We then demon- 
strate in section III that the metric and the gravitational 
slip can indeed be reconstructed without further assump- 
tions, but that the effective Newton's constant is not 
fully determined. In section IV we illustrate the practical 
consequences of this result by devising consistency tests 
for different classes of models. For general dark energy 
models in the quasi-static limit and dynamical Horndcski 
models with anisotropic stress, we are able to formulate 
general constraints that relate different scales at fixed 
redshifts and that can be tested by observations without 
the need to parametrize any free functions. For mod- 
els that have no anisotropic stress, we propose a slightly 
more complicated procedure using k- essence as a worked 
example. We are again able to test whether this model 
class is compatible with the data, and if it is, a recipe al- 
lows one to reconstruct the model's free functions with- 
out parametrization. We summarize the definitions of 
the variables used and the status of their observability in 
table I located on page 11. 



II. ASSUMPTIONS AND NOTATION 

In Ref. [18], we investigated which properties of the 
constituents of the late universe are observable without 
making prior assumptions as to the model of dark energy. 
We started off with a minimal set of assumptions, which 
we require to be able to interpret our observations at all: 

• The geometry of the universe is well described 
by scalar linear perturbations in a Friedmann- 
Lemaitre-Robertson- Walker metric with scale fac- 
tor a(t), 



ds 2 



(1 + 2*)dt 2 + a 2 {t)(l + 2<P)dx 2 



(1) 



We mostly neglect spatial curvature, but since it 
is observable, its existence would not change the 
principal results of this paper. The parameter flko 
would enter some of the equations, but given how 
small it seems to be, it would not affect them in 
any significant way. 

• The matter content is pressureless. We do not dif- 
ferentiate between dark matter and baryons. We 
neglect radiation. We assume that the equivalence 
principle holds and therefore that both these com- 
ponents follow geodesies of the same metric. This 
allows us to assume that there is no velocity bias 
between dark matter and galaxies, although we will 
comment on what breaking this assumption would 
mean. 

• The galaxy distribution is related to the matter 
distribution through a potentially time- and scale- 
dependent linear bias, <5 ga i = b(k,a)S m . We make 
no assumption about any particular form of the 
bias function. 



We have then assumed the idealized situation in which 
we can perform current and future cosmological measure- 
ments to as good an accuracy as is required for the range 
of scales and redshifts we are attempting to character- 
ize: supernova luminosities, the distribution of galaxies, 
baryon acoustic oscillation measurements and weak lens- 
ing shear. 

We have argued that without an a priori assumption 
of a particular model or parametrization for dark energy, 
there is no way of knowing the evolution or even the ini- 
tial conditions for the dark energy distribution. Indeed, 
one of the aims of late-universe cosmology is exactly to 
measure this, both on a background and perturbation 
level. 

We shall use the standard notation to describe the evo- 
lution of the cosmological background 



e{z) 2 = n m0 (i + zf + n M (i + zf + n 



x, (2) 

with E(z) = H(z)/Hq, the dimensionless Hubble param- 
eter, f2 m0 the contribution of matter to the total density, 
fix the contribution of dark energy and ttko the spatial 
curvature. 

For the purpose of perturbations, we can describe the 
solution realized by the Universe using a number of ra- 
tios of the gravitational potentials in the metric — $ and 
^ — and the dark matter configuration — the density con- 
trast 8 m and its peculiar velocity divergence 9 m . In all 
of the following, we denote with perturbation variables 
(e.g. $, ^, 5m, etc.) the square root of their power spec- 
trum, as is common in the literature. These are therefore 
positive-definite quantities and their ratio is well defined. 
We will always be working in Fourier space and primes 
will denote a derivative with respect to In a. 

We define the slip parameter 77, the effective dimension- 
less Newton's constant Y, and the dark-matter growth 
rate /, which are all in principle arbitrary functions of 
scale and redshift 



r)(k,z) = Y(k,z) = - 



2k 2 ^ 

30 m $m 



(3) 



f(k,z) 



H5- 



m 



where the scale k has been expressed in units of the cos- 
mological horizon, k = k com /aH and the prime denotes 
a derivative with respect to In a. Our definition for / is 
a little unusual, but it reduces to the standard one suffi- 
ciently subhorizon, whenever the gravitational potentials 
are slowly varying on the scales in question, since by the 
continuity equation, 



SL + = -3$' . 



(4) 



We can also think of Y as describing the size of the 
dark energy perturbation relative to the one in matter, 
although adjusted for anisotropic stress; at subhorizon 
scales and assuming no correlation between matter and 
dark energy fluctuations, 

7]Y 1 + cTT~ • ( 5 ) 
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One can form other ratios of all these variables, which 
are not independent. In particular, a ratio useful for us 
will be 



w 



f 



3(i + zfe m 

2HE 2 k 2 ^ 



(6) 



One should always bear in mind that these ratios in prin- 
ciple describe the particular solution realized by our Uni- 
verse rather than necessarily representing fundamental 
properties of the dark energy. Both Y and n are unity 
on sub-horizon scales in ACDM, while / ~ f^ 55 m the 
concordance model. All of these variables are scale inde- 
pendent subhorizon in ACDM. This is not the case for 
most dark energy models. 

It is important to stress that even measuring all 
these variables does not give all the information on the 
configuration of dark energy. In principle, there could 
be many internal degrees of freedom which do not have 
an effect on the gravitational field produced (essentially 
isocurvature modes). 

In Ref. [18], our focus was on determining the proper- 
ties of dark energy by attempting to extract from obser- 
vations its impact on dark matter. We first recovered the 
result of Ref. [19], demonstrating that no combination of 
geometrical probes of background evolution can ever split 
the contribution to total energy density of the universe 
of matter and dark energy. These probes can at best 
recover the Hubble parameter as a function of rcdshift 
up to an overall normalization, E(z), and independently 
measure any spatial curvature flko- Such properties as 
the dark energy equation-of-state parameter, w, or f2 m o 
are not observable, but can only be obtained by making 
assumptions on the form of w or performing a particular 
paramctrization. This degeneracy has been dubbed the 
dark degeneracy. Measurements of the local expansion 
rate Ho have been performed [20]. However, they are 
sensitive to the local cosmic variance and therefore do 
not necessarily reflect the averaged expansion rate today 
[21]. 

We then extended this study to perturbations. We 
have demonstrated that weak lensing and galaxy distri- 
bution observations can be reduced to a measurement of 
three variables as function of scale and rcdshift 



A = GbS m fi , R = GfS mfi 
L = Q Ia0 GY(l + r,)5 mfi , 



(7) 



where G is the growth function normalized to one to- 
day, 5 m fi is the scale-dependent density contrast of dark 
matter. In Ref. [18] we used the notation <5 mj o = cs^t.o 
to make explicit the fact that there is an unknown dis- 
tribution of dark energy perturbations present which af- 
fects gravitational measurements. As we have previously 
argued, the amplitude of the density contrast today de- 
pends on the whole history of the common evolution of 
the unknown dark energy and dark matter, and the initial 
conditions in both these components. It is thus impossi- 
ble to measure S m fl without a knowledge of the bias b: the 



scale-dependence of S m fi is unknown without a knowledge 
of the dark energy model and there are no measurements 
of the amplitude cr 8 that would be independent of the DE 
model (see the discussion in section V). Thus we have 
argued that the only model-independent observables are 
ratios of the measurements (7): 



P 1 =p = R/A = f/b, 

P 2 =E G = L/R = Q m0 Y(l + rj)/f . 

P 3 = R'/R = f + f'/f , 



(8) 



where for Pi and Pi we have provided also the notation 
more frequently used by the community [22]. It is pos- 
sible to measure each of these observables as a function 
of scale and redshift. Any other such model-independent 
observable can be re-expressed as a combination of the 
above and their derivatives. Indeed, a combination that 
appears frequently is 



T 



p?. 



p' 

Pi 1 



(9) 



which is not independent of the others. We will refer 
to the observables Pi as the primary model-independent 
observables. The observable P± contains the bias which 
cannot be measured in a model-independent manner and 
therefore we will not use it any further in our discussion 
(see the discussion of measurability of the bias in sec- 
tion V). We should stress that these results, for example, 
imply that the growth rate of dark matter / cannot be 
measured without picking a particular model or at least 
a paramctrization for dark energy. 

We will now approach the question of observability 
slightly differently than above, in a way that may be 
somewhat more physical. We will show that what is ac- 
tually observable in a model-independent manner are the 
gravitational potentials, i.e. the metric. What is not pos- 
sible to do using the observations we can perform is to 
split the sources of the curvature into a dark matter and 
dark energy part. 



III. DIRECTLY OBSERVABLE PROPERTIES 

We are going to begin this section by demonstrating 
that using weak lensing (WL) and redshift-space distor- 
tions (RSD) we can in principle reconstruct the metric 
in which the galaxies and light propagate as a function 
of redshift and scale. 

Weak lensing is a direct probe of the lensing potential 



3(1 + z)< 
2E 2 



L. 



(10) 



Thus a tomographic weak lensing map should really be 
thought of as a map of the lensing potential as a func- 
tion of rcdshift and scale [23] . It is not actually a direct 
measurement of the underlying distribution of dark mat- 
ter, which can only be extracted by accounting for the 
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unknown perturbations of dark energy also contributing 
to the potential. 

The above statement is quite natural and unsurprising. 
We will now show that RSD can be thought of in the 
same way. Traditionally one maps the velocities of the 
galaxy field onto the underlying dark-matter distribution 
through the bias and relates the RSD to the dark matter 
growth rate. 1 This is not necessary, however. Funda- 
mentally, the Kaiser formula [24] for the redshift-space 
galaxy density <5| al is a statement about a correction to 
the real-space galaxy number density <5 ga i resulting from 
the peculiar velocities of the galaxies ga \ [25], 



Vi(fc,z) 
H 



(11) 



with fi the direction cosine. Now, rather than relating 
this velocity to dark matter evolution, we can make a 
choice that allows us to avoid any assumption about the 
velocity bias: the galaxies are test particles propagating 
on geodesies of the metric, whatever the actual source of 
the curvature, 



Sal 



(12) 



with k = k com /aH. We can now integrate, neglecting 
the integration constant, to discover that from the an- 
gular dependence of the two-dimensional galaxy power 
spectrum we can extract, 



A = 5. 



gal 



R= — 



H 



= - (a 2 H) 1 j a 2 Hk' 



*d In ( 



(13) 

This approach allows us to see that linear RSD are really 
a measurement of the gravitational potential "J which 
accelerates the galaxies 



E' 

fc 2 * = R' + R ( 2 + — 
hj 



(14) 



Since all the terms on the r.h.s. of Eq. (14) are ob- 
servables, we can use RSD to measure the gravitational 
potential as a function of redshift and scale. Com- 
bining this with the information from the lensing (10), 
allows us in principle to map out the metric on our 
past light cone. This procedure will break down on 
scales where the galaxies are no longer propagating on 
geodesies, i.e. where interactions become significant or 
non-linearities of the metric invalidate the linear approx- 
imation. In principle, one should extend this simple esti- 
mate for relativistic effects which will grow in importance 
closer to the horizon, as in Refs [26-30] . 

In models where both the baryons and the dark 
matter feel the same forces (in particular, not in 
e.g. coupled quintessence [31]), we can assume that there 



1 See e.g. [24, 25]. This is how we proceeded in Rcf. [18]. 



is no velocity bias between galaxies and matter, thus 
$gai = #m- Through the dark-matter continuity equation 
(4), we can recover the standard result (7) that R = 5' m 
for sufficiently slowly varying potentials, which was the 
basis of our discussion in [18]. When this assumption 
is relaxed, which implies a violation of the equivalence 
principle, 9 m /Ho is no longer an observable because 
of the a priori unknown velocity bias. It will then no 
longer be true that RSD directly give the derivative of 
the dark matter density contrast, S' m . 

We will now attempt to characterize the dark energy, 
which will require us to return to the split between the 
dark matter and the dark energy. Indeed, the total infor- 
mation that could potentially be extracted from gravity- 
based experiments on the dark energy configuration real- 
ized in our Universe is encoded in the ratios (3). We will 
now show that there is a fundamental difference in our 
ability to extract the slip parameter r\ from observations 
and the effective Newton's constant, Y or the growth rate 
/• 

Firstly, since through Eqs (10) and (14) we can obtain 
the potentials, we can directly measure the slip parame- 
ter by taking the appropriate ratio 



r) = = 



3P 2 (l + z) a 



2£ 2 (P 3 + 2 + f ) 



1 



(15) 



This is the same result as we obtained in Ref. [18]. How- 
ever, in the current approach, we have not had to make 
any assumptions on the behavior of dark matter, or even 
its existence. Given sufficiently precise measurements of 
the galaxy distribution and weak lensing covering com- 
mon redshifts and scales, the slip parameter can be mea- 
sured redshift by redshift, scale by scale, without having 
to depend on a better- or worse-motivated parametriza- 
tion and constraints that are obtained through a total 
integrated signal. 

Another result that will be useful in what follows is 
that, since the potentials are observable, their evolution 
can also be measured. We can define an observable akin 
to the growth rate for DM, the evolution rate for the 
Newtonian potential, 



L H 

T =— = - 1 

f L l + i] 

$' n' V 

— = r + — = — 1 — 

<f> r\ L 7} (1 + rj) 



(16) 



or equivalently, the evolution of \f r can obtained purely 
from RSDs in terms of P3 by differentiating Eq. (14). 
Moreover, it is also possible to measure the variable w 
defined in Eq. (6) 



1 



Pi 



(17) 



Again, the observability of w essentially depends on 
there being no equivalence principle violation, since 
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otherwise we do not have access to 9 m . 

The measurement of anisotropic stress is a very impor- 
tant discriminator between various classes of dark energy 
models. However, in order to characterize the configura- 
tion of dark energy in our universe, we also have to map 
the effective Newton's constant, Y, which describes the 
relative size of the energy-density perturbations of dark 
energy and matter, Eq. (5). Unfortunately, it is impossi- 
ble to construct a combination of observables (8) which 
would allow one to obtain a model-independent measure- 
ment of Y. This is a result of the fact that there exists 
no cosmological probe which would be able to determine 
the dark matter density <5 m without a prior assumption 
of a particular dark energy model. As we have argued 
above, both lensing and RSDs probe the metric, with 
RSDs at best providing a measurement of the derivative 
6'^, provided the potentials vary sufficiently slowly and 
the equivalence principle be not violated. The best that 
we can do is to differentiate the definition of Y, Eq. (3), 
to obtain the differential equation 

y + n BlQ mY = i + r . (is) 

where we have dropped a term related to irrelevant 
subhorizon. The coefficients in Eq. (18) are observables 
(assuming the equivalent principle), but we have only 
succeeded in obtaining a first-order differential equation 
for Q m oY. Given measurements of observables (8), equa- 
tion (18) can be integrated. However, an initial condition 
must be supplied, effectively the value of 5 m at the initial 
time, which is not an observable. 

The above limitation is again a result of the dark de- 
generacy, here appearing on the level of perturbations. 
There is no gravity-based experiment which could tell us 
in what part the local gravitational potential is sourced 
by matter and in what part by dark energy. In princi- 
ple, one could hope to measure the mass density of dark 
matter particles by using information from direct and 
indirect detection experiments. However, since no ex- 
periment can currently determine the physical density of 
dark matter on the scales of interest, the split between 
dark matter and dark energy is fundamentally arbitrary. 
We can always redefine the dark energy model to include 
a part of the dark matter. 2 

A part of the dark degeneracy is related to our lack of 
knowledge of f2 m o and implies that Q m o can be absorbed 
into the definition of Y through Y — > tt m oY. However, 
this is not sufficient to fully specify the split between dark 
matter and dark energy: we still have the freedom to sup- 
ply an arbitrary initial condition for (18). For example, 



2 See for example Ref. [32] where a model with no dark matter, 
just dark energy, was proposed. The dark energy has sound speed 
c 2 = and can, for example, replicate exactly the phenomenol- 
ogy of the background and dark-matter linear perturbations in 
ACDM. 



one could imagine a situation where at some particular 
redshift there are no dark energy perturbations but only 
those for matter, i.e. Y = 1. Once this choice has been 
made, the observations in combination with Eq. (18) tell 
us how Y evolves as a function of redshift. However, any 
initial choice is as good as any other. Eq. (18) essentially 
evolves the values of Y to the next timeslice as a result of 
our assumption that the part of the universe that we have 
designated to be dark matter is prcssurclcss and follows 
geodesies of a metric we have mapped. 

As a result of this degeneracy, the best that we can do 
is to propose a behavior for Y, resulting from some sort 
of model with some free parameters, and then ask the 
question whether the observations are consistent with 
the evolution mandated by Eq. (18). Eq. (18) should 
therefore be thought of as a consistency relation for Y 
which can be used to rule out particular models but not 
to measure this quantity. We will give an example of 
such a procedure in section IV. 

The reason for this fundamental limitation is our lack 
of knowledge of the density bias. The measurement 
of the amplitude of the galaxy power spectrum only 
describes galaxy number counts; given a measurement 
of the bias for some particular tracer, the dark matter 
density contrast could also be reconstructed. Methods 
for the measurement of bias have been proposed and 
used, but they assume various properties of ACDM in 
a fundamental way. They cannot be used to constrain 
dark energy models departing from the cosmological 
constant in general (see section V for a discussion). 

In light of the above results, we can reinterpret 
slightly the measurements that we are performing 
in ACDM. The weak-lensing and galaxy-distribution 
observations are in fact sensitive to the perturbed 
geometry, since they involve tracers that propagate on 
geodesies. It is only because the cosmological constant 
is not perturbed that these potentials can be ascribed 
directly to dark matter and the standard language 
recovered. It should be intuitively clear that without 
the knowledge of dark energy perturbations we must 
have arrived at the model-independent result we have 
described: only geometrical and not dynamical proper- 
ties can be measured by tracers propagating on geodesies. 

The overall conclusion is that without particular as- 
sumptions on the model of dark energy, cosmological 
observations allow at most for the mapping out of the 
dimensionlcss Hubble parameter E{z) as a function of 
redshift and the two gravitational potentials as a func- 
tion of both redshift and scale. No information about the 
relative distribution of dark energy and dark matter can 
be obtained without an assumption of a particular model 
or at least some paramctrization. 
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IV. TESTING SCALAR-FIELD MODELS 

In this section, we will discuss our ability to use the 
observable properties described in section III to constrain 
or even exclude generalized scalar-tensor models of dark 
energy. 



A. Quasi-Static Dark Energy 

As we have already discussed in Ref. [18], perturba- 
tions in scalar-tensor dark energy models frequently can 
be approximated to evolve in the quasi-static approxi- 
mation, where the dark energy follows the dark matter 
perturbations and the time derivatives are negligible (see 
also Ref. [33]). For the Horndeski Lagrangian [34, 35], 
which is the most general scalar-tensor theory involving 
no more than second derivatives, this limit was derived 
in Ref. [36] . This derivation assumes that the only scales 
relevant to the problem arc the Jeans length (determined 
by the sound speed) and the Compton wavelength, deter- 
mined by the effective mass of the scalar perturbations. 
Dispersion relations arc in principle more complicated 
in these models, but apart from the mass, the sound 
speed does provide the smallest scale by the virtue of 
being defined in the k — > oo limit. Under the assumption 
that the speed of sound is close to that of light and the 
scales larger than the Jeans length lie outside of the scales 
probed by the observations, the effective Newton's con- 
stant and the slip parameter (3) in these general scalar- 
tensor models take the form 

Y = hl {TTkth;)> ^ = HiT^J' (19) 

where the functions hi are purely functions of rcdshift 
determined by the Lagrangian describing the dark energy 
model (sec Ref. [18] for the explicit form of his). 

The quasi-static limit is really the requirement that the 
dark energy perturbations follow the dark matter ones 
in a very constrained way. Dark matter evolves on the 
geodesies of the combined gravitational potential while 
the dark energy perturbation must follow the very precise 
prescription defined by Eqs (19), without any dynamics 
of its own. 

One can use Eq. (15) to test the consistency of the ob- 
served universe with the slip parameter r\ given by the 
form (19). As we have described in Ref. [18], given mea- 
surements of the observable (15) at more than three dif- 
ferent scales per redshift, we can test whether the data 
are consistent with a description for 77 of the form (19). If 
not, then the anisotropic stress cannot be described as a 
manifestation of a scalar-tensor theory in the quasi-static 
limit. 

Since dark matter is supposed to move in a known way 
on geodesies of the metric and we can map this metric 
out, the consistency of the form of Y given in the quasi- 
static limit can be tested, despite the non-observability 



described in section III. Inserting the form (19) into the 
consistency relation (18) allows one to ask, redshift by 
redshift, whether the quasi-static form for Y is consistent 
with the observational data, i.e. whether the simple scale 
dependence for Y driven by the coefficients hi is suffi- 
cient to explain the observations. As a result of taking 
the time derivative in (18), the number of free param- 
eters increases to seven (including £l m o) , thus a larger 
number of measurements at each redshift is required to 
overconstrain the system, but it is in principle not any 
more difficult than for the anisotropic stress. 

The main takeaway from this discussion is that if the 
quasi-static limit can be assumed, the dark energy con- 
figuration can be tested by looking for a particular scale 
dependence of observables redshift bin by redshift bin, 
rather than as a single integrated fit to the observations 
with a particular parametrization chosen for the func- 
tions hi. 



B. Beyond the Quasi-Static Limit 

The quasi-static limit considered in section IV A is 
somewhat unsatisfactory, given that there are not nec- 
essarily a priori reasons why the DE should evolve so as 
to follow the DM perturbations so closely. Y specifics 
the size of the dark energy density perturbation relative 
to the dark matter one. It should therefore be intuitively 
clear that on any single spatial hypcrsurface, one can in 
principle arrange the dark energy density in an arbitrary 
manner and therefore there can be no restriction on the 
form of Y if the dark energy can be perturbed freely. 
By choosing the quasi-static limit, we give up this free- 
dom, since all the relative velocities between DM and 
DE must vanish and the distribution of the DE is purely 
determined by the density perturbation of dark matter. 

Nonetheless, once the initial conditions are specified, 
Y will have to evolve in a manner consistent with the 
equation of motion for the scalar or, equivalcntly, the 
conservation of the DE energy-momentum tensor, what- 
ever the model behind it. Since the metric is mapped 
out, and the dark matter propagates on its geodesies, it 
is possible to test whether the evolution of the dark en- 
ergy part is consistent with a particular class of models. 

Another requirement arises in situations in which non- 
zero anisotropic stress is detected. Its measurement 
through Eq. (15) gives additional information on the dis- 
tribution of the scalar-field perturbations. This then has 
to be consistent with the distribution implied by the ef- 
fective Newton's constant Y on that spatial slice. The- 
ories with anisotropic stress are in general more compli- 
cated than those without and contain more parameters. 
However, in this way we can have an additional handle 
with which to constrain them. 

We will now discuss a couple of examples to illustrate 
how one might go about constraining models using the 
model-independent variables. 
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1. Consistency with k-essence 

Here, we will turn to the simplest class of scalar field 
models, those described by the Lagrangian £^ = K(cj>, X) 
[37-39]. The energy momentum tensor possessed by this 
kind of dark energy has perfect-fluid form and at linear 
level all the properties are described by the equation of 
state w and the sound speed c s [40]. This example will 
provide for sufficient complexity in order to demonstrate 
the logic of our method and the fundamental limitations 
of constraining the dark energy model space as a result of 
the non-observability of Y, while having the advantage of 
being familiar to a wide audience. Note that all uncou- 
pled quintessence models are contained within this class 
as are perfect- fluid models, provided only scalar pertur- 
bations be considered. 

We shall not give the full Einstein equations here, but 
it will suffice to say that the combined k-essence and dark 
matter energy-momentum tensor (EMT) depends on the 
perturbation variables schematically as 

ST° D S<j>, 54> , S m , (20) 

sifDS<f>,e m , 

STj - l /35ftTg = 0, 

with 5(f> the perturbation of the k-essence scalar and 
5m, #m describing the DM density contrast and the di- 
vergence of the peculiar velocity. This model does not 
allow for any anisotropic stress and therefore it could be 
immediately excluded if, even at just one redshift and 
scale, T) ^ 1. We will assume that no such detection was 
made and therefore we will take $ = — We can then 
use the combination of the Hamiltonian and momentum 
constraints [i.e. the (00) and (Qi) Einstein equations], to 
eliminate the scalar-field perturbations S<p and 8(j> in the 
(ii) Einstein equation, obtaining the exact equation for 
the evolution of the gravitational potential in k-essence 
models in the presence of dark matter 

*" + (4 + § + ZcA ^+U + 2^ + 3tA (21) 



c 2 fc 2 * 



-Q m (c 2 6 m + 3(c 2 - cDk^H-Hm) 



where the adiabatic sound speed, 

2 = Px 6E'/e + 2(e'/e) 2 + 2E"/e 
a ~ Px ~ 9O m0 £- 2 (l + z) 3 + 6£y B ' 1 ; 

is fully determined by the observable expansion history 
up to n m o • 0-m can be neglected sufficiently deep sub- 
horizon, although that term it is actually observable as 
0m/H = —R and can also be included in the calculation. 
Equation (21) is quite standard, see the closely related 
result in e.g. [41, Eq. 7.51] 

Apart from the parameters of this model — the constant 
flmo and the sound speed c s , which can be a function of 



time — all the quantities on the left-hand-side of Eq. (21) 
are observable. What is not observable is d m , as we have 
explained in section III. Since, at any one redshift slice, 
the dark-matter configuration can in principle be arbi- 
trary, Eq. (21) does not by itself provide a constraint on 
the theory space. However, we can think of Eq. (21) as 
a measurement of Y given the assumption that the dark 
energy model belongs to the k-essence class: 

r 2 k 2 

Y = ^- 

[F + (r + 3 + E'/ E + 3d) (1 + r) + E'/ E + c 2 k 2 ' 

(23) 

where we have used Eq. (16) to replace derivatives of "J 
with the observable T. So given observations of T and 
E, Y [c^ , H m o] is a functional on that data depending on 
the parameters f2 mu and c 2 which outputs a function of 
scale and time. 

However, we should stress that we can always find some 
Y given any data and given any choice of parameters. 
We have therefore thus far only obtained a model for Y 
in the spirit of that provided by (19). It is a little more 
complicated, in that it depends not only on parameters 
equivalent to the h^s, but also on observable data. If 
the observable T shows no scale dependence, in principle 
we have a very simple model for Y, not dissimilar to the 
quasi-static one. If on the other hand, the observations 
do show scale dependence, Y could be a very complicated 
function of scale. 

We must now test whether this model for Y is consis- 
tent with the observations using the consistency relation 
(18), just as we did in the quasi-static case. This means 
that k-essence can only be a good description of the dark 
energy if there is no anisotropic stress observed and also 



Y 



2i7inoF 

P 2 



1 



(24) 



at every redshift and every scale. This consistency rela- 
tion must be valid given just one global parameter f2 m o 
and at each redshift the sound speed c 2 and its derivative. 
Just as in the case of the quasi-static limit, at any one 
redshift, given measurements at four scales it is in prin- 
ciple enough to exclude such a model for Y and therefore 
k-essence as the mechanism for dark energy, assuming 
the existence of a perfect data set, as we already de- 
clared. 3 Note that our approach is independent of the 
initial conditions. 

If given the assumption of k-essence, the consistency 
relation is satisfied at all redshifts, then the above 
measurements can be used to determine the value of the 



3 In fact, this would exclude any single perfect fluid as dark en- 
ergy, since k-essence is equivalent to perfect-fluid hydrodynam- 
ics. Multiple perfect fluids or fluids with internal degrees of free- 
dom have in general more complicated pressure perturbations, 
which would naively appear as more complicated scale depen- 
dence (e.g. [42]). 
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sound speed at each redshift observed. This allows a 
non-paramctric determination of this physical property 
that in general is an arbitrary function of time. 

It is worth mentioning the limit c 2 = [32, 43] , which is 
not quasi-static on any scale. In such a case, the dust and 
dark energy perturbations become indistinguishable and 
the entropy perturbation source in Eq. (21) disappears, 
yielding 



5 + 35 » 



3 + 2^ + 3ci 



(25) 

where we have redefined the adiabatic sound speed to be 
that corresponding to the total EMT, c 2 = ftx/{f>x + 
p m ), i.e. it is now a function purely of the background 
geometry. In this limit, there is no dependence on the 
unobservable S m and all scale dependence disappears. A 
measurement of Y is not possible in a fundamental sense, 
since there is no difference in the properties of the DE and 
the DM perturbations. The dark degeneracy is complete. 
The physics of the linear perturbations is effectively that 
of a single generalized dust collapsing on a background 
with some equation of state [32]. However, for such a 
model to be a valid description of the observations, the 
data must satisfy the constraint 

r'+ (r + 3 + 2^ + 3cl\ (r + i) + ^ = o (26) 

at every redshift and every scale with no free parameters. 
This constraint reduces to exactly that of ACDM when 
the equation of state for the dark energy becomes — 1. 

Models with kinetic gravity braiding (KGB) [44, 45] 
are the most general class of scalar-tensor theories with 
a single scalar which do not have a direct coupling to 
gravity and therefore do not have anisotropic stress de- 
spite being imperfect fluids [46]. The KGB equivalent of 
Eq. (21) would feature more scales, however the prescrip- 
tion for constraining this class of models would not differ 
from the k- essence example presented here. 



2. Consistency with general Horndeski theories featuring 
anisotropic stress 

The most general second-order scalar-tensor theory is 
described by the Horndeski Lagrangian [34, 35]. The 
structure of the energy-momentum tensor for these the- 
ories is much more complex than for a model such as 
Y-essence. In particular, the Lagrangian features a non- 
minimal coupling to gravity. This means that the EMTs 
for any such models will feature, in general, anisotropic 
stress related to the perturbations of the scalar. 

Horndeski models in the presence of dark matter have 
EMTs which have the following dependence on the per- 
turbed fields when only linear perturbations are consid- 



ered: 

5T^S(t>,S4>,S m , (27) 
5T?D5<t>,64>,6 m , 
STj - yaSjSTg D aS(f> 7 

6Tf d 5<j), 5<f>, 8(j) , 

in addition to dependence of all the components on the 
gravitational potentials $ and \& and their time deriva- 
tives. Here a is fully determined as a function of the 
Horndeski free functions. 4 As compared to k- essence, 
we have the already mentioned anisotropic stress, which 
is always proportional to 8<j). The perturbation 5(f) is 
present in the (Oi) components of the EMT, while the 
pressure perturbation depends on 8<j). In addition to the 
above, we have the equation of motion for the scalar field, 
which is an equation for 6(j> in terms of all the other vari- 
ables. The full expressions can be found in Ref. [36]. 

We can clearly see that the only way to suppress the 
anisotropic stress is to cither make the coupling a very 
weak, i.e. effectively make the Horndeski terms non- 
minimally coupled to gravity irrelevant for the dynamics 
of dark energy, or to suppress the scalar perturbations, 
which is only possible if the scalar is very massive and 
not evolving, i.e. it just contributes vacuum energy [47]. 

Using the equation of motion for the scalar, we can 
eliminate 5<p from the pressure equations. Then, using 
the Hamiltonian, momentum and anisotropy constraints 
[i.e. the (00), (Oi) and (ij) Einstein equations], we can 
eliminate three further variables: the scalar perturba- 
tions 6(j>, 6(j> and the dark-matter density 5 m . This is one 
extra variable as compared to models without anisotropic 
stress, where the anisotropy constraint equates the two 
potentials but is not dependent on the configuration of 
the dark energy degree of freedom. The remaining form 
for the (ii) Einstein equation is very simple, 



$" + + a 2 *' + (a 3 + a 4 k 2 ) $ + (28) 
+ (a s + a e k 2 ) ($ + #) = aa 7 n m k,- 2 m . 



where the parameters on are functions of time alone and 
are fully determined by the Horndeski Lagrangian. Their 
exact form will not be useful here. 

We have kept the dependence on a explicit to show 
that in the limit of vanishing anisotropic stress, a — > 0, 
Eq. (28) is not an alternative dynamical equation for the 
evolution of the potentials To obtain the evolution 
equation, one would have to eliminate ^ for $ and would 
obtain a very complicated equation with many scales em- 
bedded in its coefficients. One may worry that Eq. (28) 
provides a different version of this evolution equation 



4 It is equivalent to the parameter B-j in Ref. [36]. It is only non- 
zero when there is a direct coupling of the scalar to gravity in 
the action. 
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than that obtained for theories with no anisotropic stress. 
This is not the case: switching off the anisotropic stress 
reduces Eq. (28) to an equation which enforces $ = 
thus providing no additional information. 

Eq. (28) is exact at all scales where linear perturbation 
theory is valid. In Eq. (28), all the variables are model- 
independent observables. Using the extra information 
from the anisotropic stress constraint, we have eliminated 
<5 m , which is not observable, and thus we can obtain a 
consistency relation between pure observables required 
by all Horndcski models with anisotropic stress, 

r)T' + 7]" + T(i 1 T + 2t 1 ' + a lV + a 2 )+ (29) 
+ airf + a 3 n + a 5 + k 2 (6:477 + a 6 ) = ajw . 

where we have redefined the on to absorb factors of er, E 
or f2 m , or to combine coefficients into single variables. 

A nice simple example of the relation (28) can be ob- 
tained in the case of the f{R) class of dark energy mod- 
els, which are a subclass of the Horndeski lagrangians 
and therefore feature much less freedom [48]. Following 
the algorithm described above, one obtains as the exact 
linear result 

- + ( 4 + J) + ^ ( m c + 2fc2 ) (30) 

+ K m2c_6 ( 2+ S +fc2 ) t=0, 

where m 2 -, = J \r/2H 2 'j \rr is the Compton mass of the 
scalar degree of freedom in the units of Hubble. We must 
have m c 3> 1 in order to satisfy Solar-System constraints 
[49]. At small enough scales, k 3> mo, we can recover 
the standard quasi-static result in f(R) that 2$ + ^ = 
0. Again, we can express the above as a null test on 
observables that needs to be satisfied at all scales where 
linear theory is valid 

r}T' + jj" + r Lr + 27/ + 5 + + ( 4 + ? /+ 

m 2 , s ( E'\ k 2 , 
+ ^(^-l)+6(2+— J + — (277 - 1) = . 

(31) 

This relation has only one free parameter, m^, which 
can be freely adjusted at each redshift to fit the data. 
If this is not enough to satisfy the above test, then 
f(R) is not a general enough theory to account for the 
observations. 

The relation (29) [and (31)] is valid at all linear scales, 
without the need to involve the quasi-static limit or 
choose initial conditions. It is enough to make a detection 
of non- vanishing anisotropic stress at even just one scale 
and redshift to then require that the data satisfy (29). At 
any given redshift, the a, coefficients are just numbers, 



independent of scale. The observables, T, 77 and w are in 
principle complicated functions of scale and Eq. (29) is a 
non-linear scale-dependent function on these data. Per- 
forming observations for at least eight values of k one can 
form an overconstrained system and rule out the Horn- 
deski model, or confirm it while measuring its parameters 
ai_7, redshift by redshift without prior parametrization. 



V. DISCUSSION AND CONCLUSIONS 

Dark matter is dark. This statement is obvious and 
yet the discussion of cosmological observations is usually 
framed as a discussion of the measurements of the dark- 
matter distribution. This is of course completely natural 
in the case where there are no dark energy perturbations, 
i.e. within the framework of ACDM cosmology. However, 
the moment that we relax this assumption, as we must 
do when investigating dark energy models different from 
a cosmological constant, this connection between mea- 
surements and the dark-matter distribution is no longer 
simple, and can to all intents and purposes disappear. 

In this paper, we have stressed that the probes we hope 
to use to constrain general dark energy models depend 
on tracers which propagate on geodesies: light for weak 
lensing and galaxies for RSD. This allows us to map out 
the metric on our past lightcone (the potentials <f> and \I>) 
through which these tracers fall freely. These arc the real 
physical observables of the combination of the cosmolog- 
ical probes. The relation of these potentials to the dark- 
matter distribution is then a model-dependent statement, 
which just happens to be an identity for ACDM. 

Given this limitation, we should stress that the cos- 
mological probes still can provide an enormous amount 
of information. Since we can reconstruct the potentials, 
the slip parameter 77 of the matter/light Jordan- frame 
metric is a model-independent observable, given overlap- 
ping WL and RSD measurements. We do not need to 
parametrize it and we can answer the question of whether 
anisotropic stress is at all necessary, directly from ob- 
servations, without any further modeling. This test has 
the power to immediately eliminate very large classes of 
models, whatever its result. If 77 ^ 1 at even one redshift 
and scale, then models without non-minimal coupling to 
gravity in the baryons' Jordan frame can be thrown away. 
If 7] = 1 everywhere, then such couplings of dark energy 
to gravity could in principle still be there, but they are 
so small that they cannot at all influence the dynamics 
of the universe and therefore can be neglected. We have 
summarized the most important dynamical variables and 
their observability status in Table I. 

We have also shown that the evolution rate of the po- 
tentials can be mapped out. On the other hand, the 
fact that we cannot observe the dark matter density am- 
plitude means that the effective Newton's constant Y is 
not an observable in a model-independent setting. Given 
these model-independent observables, we should refocus 
the predictions of perturbation theory in dark energy 
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models to their impact on these gravitational potentials 
rather than on the dark matter, contrary to the dominant 
approach in the discussion today. 

Nonetheless, we have demonstrated in section IV how 
to construct tests of very general classes of models of 
dark energy, even away from the frequently employed 
quasi-static limit. The distinct advantage of this 
approach is that rather than parametrizing the free 
functions of a class of models first, then evolving the 
predictions of these models using modified codes and 
simultaneously fitting to data, our method is capable of 
testing/constraining classes of models without any prior 
paramctrization and without any assumptions on the 
initial conditions. These tests involve performing mea- 
surements of the observablcs at multiple scales, rcdshift 
by redshift (with the number of data points required 
determined purely by the number of free functions of 
time in the particular class of models). In order to 
carry out this test it is not necessary to modify codes 
for each class of models, but just to apply appropriate 
transformations to the observed data. In other words, 
our approach calls for exploring the space domain rather 
than the time domain of dark energy. 

We have mostly assumed that there is no velocity bias 
between dark matter and the galaxies. Thus a mea- 
surement of galaxy peculiar velocities 6 g& \ via RSD is 
automatically a measurement of the dark-matter pecu- 
liar velocity 9 m . Statistical origins of the velocity bias 
notwithstanding [50], there could in principle be a non- 
vanishing bias between these two velocities. Indeed, this 
happens in any model in which there is a fifth force act- 
ing on dark matter, which does not couple to baryons or 
light, i.e. models with a violation of the equivalence prin- 
ciple such as coupled quintessence [31]. This fifth force 
is a new source of acceleration for dark-matter particles, 
causing their peculiar velocities to deviate from that of 
the galaxies. The effect of this is to introduce a (scale 
and time-dependent!) velocity bias. Our measurements 
of WL and RSD still map out the same metric potentials 
for our tracers, but these are not the complete potentials 
that the dark matter feels. Allowing for such deviations, 
we completely lose the connection between the dark mat- 
ter growth rate and RSD of the galaxy power spectrum. 

Summarizing the above, we can say that RSD and 
WL map out the Jordan-frame metric for galaxies and 
light, but provide no direct information on the Jordan 
frame of the underlying dark matter. In principle, one 
could use different tracers with different baryon fractions 
(cluster and galaxy power spectra) to map out the 
Jordan frames for these two classes of objects and obtain 
some information on the differences in the gravitational 
potentials they experience (in the spirit of e.g. [51]). 
Interestingly, galaxies are not just baryons, but are 
dynamically coupled to their dark-matter halos. This 
means that the Jordan-frame metric for the galaxies is 
not necessarily the same as the one for baryons, since 
the galaxies will at least partly feel the fifth force acting 



on dark-matter. 

Another interesting implication of our approach is 
that potential non-linearities of dark energy, which can 
appear on very different scales to those in dark matter 
as a result of the generic existence of screening mecha- 
nisms, are not deadly to the measurements of the linear 
potentials <& and 'J. For tracers to move on geodesies 
determined by these potentials, their gradients must 
remain small, i.e. the total EMT perturbations must be 
linear. So even if the DE contribution is non-linear, it 
is enough for the dark-matter density to be dominant 
and its perturbation linear for the tracers to continue to 
move on geodesies described by $ and ^ . However, since 
there will no doubt be non-linear contributions to the 
potentials, both from the DM and the DE perturbations, 
this should imply that measurements of higher-order 
correlation functions should be sensitive to the extra 
non-linearity as compared to the expected ACDM result. 

Let us pre-emptively address a few potential criticisms: 
our results and approach is driven by the limitations 
in our ability to determine the galaxy bias b as well as 
the normalization og and scale dependence of the dark- 
matter perturbations <5 m . There exist many, increasingly 
better, measurements of these quantities e.g. [52-56]. 
However, they all depend in a fundamental way on as- 
sumptions necessarily true only in ACDM and therefore 
we cannot use them. 

Although non-linearities can be used to constrain the 
bias through the galaxy bispectrum [56] , this is a model- 
dependent statement since screening mechanisms present 
in most dark energy models alter the expectations. Oth- 
erwise, a measurement of bias is usually obtained as 
a result of measuring ag, which is done using cither 
cluster counts and/or weak lensing around clusters. As 
we have already pointed out, weak lensing measures 
the lensing potential rather than the DM distribution, 
i.e. includes a dependence on the DE perturbations. 
To extract information from cluster counts, their mass 
is obtained from their gas temperature — which gives 
information on the Newtonian potential \& and not 
on the DM mass — and spherical collapse or N-body 
simulations are used, both of which depend deeply on 
the model for dark energy, see e.g. Refs [45, 57, 58]. 
All of these issues can in principle be appropriately 
modeled, but this must be done on a modcl-by-modcl 
basis. Without this calculation, we cannot interpret 
the ACDM measurement of a$ as a parameter valid for 
other models of dark energy. 

The tests that we have constructed in section IV B de- 
pend on eliminating variables through the use of the con- 
straints present in the Einstein equations. Given a class 
of models for DE, one can replace dynamical variables 
with combinations of observables. Given the assump- 
tion that the DE is a single degree of freedom, that the 
equivalence principle be satisfied and a detection of non- 
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Class 


Variable 


Key Relation 


Obs.? 


Eq. 


Comment 


Measurements 


A 


= 5 g al = 65m 


/ 


(7) 


Bias not measurable absent DE model 




R 




/ 


(7) 


RSD measure galaxy velocity 




L 




/ 


(7) 


WL shear probes lensing potential 




E 


= H/H 


/ 




Ho, fi m o not observable without DE model a 


Primary 


Pi=P 


= R/A = f/b 


/ 


(8) 




observables 


Pi = E G 


= L/R= (l + r))/m 


/ 


(8) 






Pz 


— R'/R = f + f'/f 


/ 


(8) 


Only this function of / is observable 


Physical 




R' + R(2 + E'/ E ) 


/ 


(14) 


Extract from RSD 


variables 


<1> 


3(1+z)3 T 1 >T< 


/ 




Extract from WL tomography 




S m 


5m — <5 g al/6 






Unknown without knowing bias 




O m / Ho 


Om = #gal 


/ 




Observable only given equivalence principle 


DE 


V 


EE -$/* 


/ 


(3) 


n 7^ 1 => coupling DE/gravity 


Configuration 


r 


= #'/# 


/ 


(16) 


Observable prediction of DE models 


Variables 


Y 


= -2fc 2 */3fi m 5m 




(3) 


<5 m unknown; must satisfy relation (18) 




f 


EE —Om/HSm PS 5'm/S-m 




(3) 


8 m unknown 




w 


= f/MmoY 


/ 


(6) 


Relates DM velocity to potential 



a See Ref. [18] for a detailed discussion of the observability of the background. 

Table I. Summary of variables used in this paper. Those variables marked with a checkmark as observable can be measured as 
a function of redshift and scale without an assumption of a particular DE model. Measurement of Sm requires the knowledge of 
bias, which cannot be measured or modeled without the knowledge of the DE model. The variables such as Y and / typically 
used to describe the DE configuration are therefore not observable, but predictions can be reformulated to maximally exploit 
those variables which are observable (see section IV B). 

vanishing anisotropic stress, we have enough information 
to eliminate all non-observable quantities from the evo- 
lution equations and thus are able to form a null test for 
the most general class of scalar-tensor models directly on 
the observable data. 

When there is no anisotropic stress then no such com- 
plete constraint can be formed: the dependence on the 
unobservable S m remains. In that case, the best that 
can be done is to obtain a measurement of the effective 
Newton's constant Y on the assumption that a particu- 
lar class of models of dark energy describes the Universe 
and then use a consistency relation for Y to determine 
whether the assumption was consistent with the data. 
This is a somewhat more complicated exercise, which 
we have demonstrated for perfect-fluid k- essence models, 
but would also apply to imperfect-fluid models featuring 
kinetic gravity braiding [35, 45, 46]. 

Tests of the type we have shown can in principle be 
constructed for any other class of DE models. Indeed, 
we would like to argue that one should think of the 
parametrized or effective approaches, such as those of 
Refs [13, 15-17], as providing the dynamics of the dark 
energy which can be rewritten in terms of evolution equa- 



tions for the potentials. Provided they are written in 
terms of parameters that on cosmological solutions are 
only functions of time, a test such as those presented in 
section IV would allow for putting constraints on these 
parameters, or indeed would exclude setups which posses 
insufficient operators to describe the data fully. 

An alternative way of thinking about these tests 
is that, using the measurement of <f> and 'J', we are 
essentially reconstructing the components of an effective 
combined EMT for DM and DE [59]. Now, given a 
class of models, we can extract relations between the 
configuration of the degrees of freedom and the fluid 
variables, and therefore between the fluid variables 
themselves, i.e. anisotropic stress, pressure, energy 
density. In general this is difficult, but in some classes 
of models it can be done, e.g. [12, 14]. However, it 
should be intuitively clear that if the dark energy has 
more degrees of freedom, only the adiabatic modes will 
influence the potentials. Internal (isocurvature) modes, 
since they do not affect the gravitational field of the DE, 
would only be constrainable through their impact on the 
time evolution of the dark energy. 
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We would like to stress that our method depends on 
taking derivatives of the data. Since for the purposes of 
this paper we are working in the idealized case of suffi- 
ciently good data, we have not addressed here the feasi- 
bility of this procedure in a realistic situation. We leave 
this for future work. 

We have also used a number of simplifying assumptions 
for the relation of measurements to the potentials. For 
example, the Kaiser formula (11), assumes not only lin- 
earity, but also a flat sky and ignores near-horizon effects. 
As shown by e.g. [29] , the correction involve contributions 
from, for example, weak lensing, which would pollute the 
determination of the peculiar velocities of galaxies from 
RSDs. Indeed there are also near horizon corrections to 
the weak lensing, e.g. [26]. Since our aim is to exploit 
large-scale data fully, taking into account such correc- 



tions is a natural extension of this work. 
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